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some simple analytical approximations 


Alexander A. Kokhanovsky 
University of Bremen 
Institute of Environmental Physics 
FB1, P.O. Box 330440 
D-28334 Bremen, Germany 
and 
B. |. Stepanov Institute of Physics 
70 Skarina Avenue 
Minsk 220072, Belarus 
E-mail: alexk@iup.physik.uni-bremen.de 


Abstract. This work is devoted to the derivation of simple approximate 
equations, which can be used for studies of light reflection from various 
semiinfinite weakly absorbing turbid media. It is assumed that the prob- 
ability of photon absorption in a single scattering event is small. The 
method used is largely based on the approximate summation of the 
MacLauren series for total light reflectance and other radiative transfer 
characteristics as functions of the probability of photon absorption. The 
main results comprise convenient analytical expressions for total light 
reflectance, plane albedo, reflection function, and degree of polarization 
of reflected light. © 2003 Society of Photo-Optical Instrumentation Engineers. 
[DOI: 10.1117/1.1578086] 
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1 Introduction 


The core of the radiative transfer theory is the vector inte- 
grodifferential radiative transfer equation,’ which describes 
the intensity and polarization of diffused waves in scatter- 
ing media. Generally speaking, most of polarization and 
transfer characteristics of various media can be derived us- 
ing this equation. However, the equation is not very conve- 
nient for analytical studies. Thus, there is a need for the 
derivation of simple analytical equations, which can de- 
scribe the transport of light in random media. 

A number of analytical approximations for various light 
field characteristics are discussed by van de Hulst in his 
fundamental monograph on the subject of radiative 
transfer.” They are largely based on various transformations 
and approximations of the radiative transfer equation. 

It is of interest, however, to derive analytical equations, 
well known in radiative transfer theory, but actually not 
using the radiative transfer equation explicitly. These ap- 
proximate equations can be used in numerous applications 
of the radiative transfer theory, especially as far as optical 
engineering is concerned. It could be of interest also from a 
didactic point of view. 

We consider here only one selected problem of the ra- 
diative transfer theory, namely the reflection of light from 
semiinfinite media (e.g., snow fields, paints, whitecaps, 
etc.) However, other more complex cases, including finite 
scattering layers, can be handled using the same chain of 
arguments. 

The reflection of light from a semiinfinite turbid medium 
is a standard problem of radiative transfer.’ Actually optical 
properties of finite but optically thick slabs can easily be 
found if the reflection function of a semiinfinite layer with 
the same local optical characteristics is known.” This is 
why that problem was studied in great detail both theoreti- 
cally and experimentally. Especially fruitful results as far as 
reflectance spectroscopy is concerned were obtained by 
Rozenberg.*~° They were updated and brought to corre- 
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spondence with the exact asymptotic radiative transfer 
theory by Zege, Ivanov, and Katsev.’ 

Our primary task is to present a new method for the 
derivation of analytical equations for total reflectance and 
other reflection characteristics of semiinfinite turbid media. 
Note that the method proposed can be used well beyond the 
narrow scope considered here. 


2 Total Reflectance of a Semiinfinite Turbid 
Layer 

Let us start from the derivation of an approximate formula 

for the total reflectance r, of a semiinfinite layer. The value 

of r, is defined as follows!: 


1 2m 


mi2 
rs 57 i ay fi sin 2 ddd 


mi2 
x | sin2 Bod PoR (Do, O, W), 
0 


where R(®g, O, y) is the reflection function” of a plane- 
parallel turbid layer, Oo is the incidence angle, 0 is the 
observation angle, and W is the relative azimuth between 
incident and reflected light beams. The reflection function 
R(0_9,0,) is defined as the ratio of the intensity of re- 
flected light for a given scattering medium to that of an 
absolutely white Lambertian surface. By definition, we 
have for such a surface: R=1 and, therefore, r,=1. 

First of all, we note that the value of r,(wọ) for semiin- 
finite turbid media can be presented in the following gen- 
eral form: 


r(ao)= >) a09, (1) 


n=1 
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where a, are unknown coefficients and wọ is the probabil- 
ity of photon survival in a single light scattering event (so- 
called single scattering albedo'). We assume that wo is 
close to 1, which is the case for many substances in visible 
and near-infrared regions of the electromagnetic spectrum." 
Then it is more convenient to use the expansion: 


co 


r,(B)= > a„(1- B)", 2) 


n= 


where 6=1— wọ is the probability of photon absorption, 
which is assumed to be a small number. Clearly, we have 
from Eq. (2) at B=0: 


r(0)= X a. (3) 


n=1 
However, it also follows: 
r,(0)=1, (4) 


which is due to the energy conservation law. Thus, one 
obtains that 


> a,=1, (5) 


n=1 


and numbers a,, can be interpreted in terms of the probabil- 
ity theory. In particular, the value of a, gives us the prob- 
ability that the photon will be singly scattered before escap- 
ing the turbid medium. The probabilities of scattering 
events 4], 42, a3,... do not depend on each other. The 
theorem of adding independent probabilities brings us to 
Eq. (5) as well. Let us substitute the following exact expan- 
sion in Eq. (2): 


(1-p)"=> |") 6) 
j=0 J 
where 
") _ n! J 
HT Fae Á 
Then it follows from Eqs. (2) and (6): 
r{B)= È and, ry") (8) 
n=1 j=0 J 
or in the explicit form: 
œ 2 _ 
r(B)= ay ppi EZEN 


Bn(n-1)(n—2) 
= oe. 


Equation (9) can be rewritten in the following form: 
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(10) 


where we used the normalization condition in Eq. (5) and 
defined the following averages: 


n=1 n=1 


(n(n—1)(n—2))= >) n(n—1)(n—2)ay, (11) 
n=1 


and so on. Here (n) is the average number of scattering 
events in the medium. 

Equation (10) is an exact formula. Now we make some 
assumptions to simplify Eq. (10). First, we assume that (n) 
is large and, consequently, (n(n—1))~(n7), (n(n—1)(n 
—2))~(n3) and so on, are also large. Clearly, such an 
approximation is valid for weakly absorbing media (as 8 
—0). Note that for nonabsorbing semiinfinite media (n) 
—, This gives us, instead of Eq. (8): 


2 3 
r= 1 pln) + (n) dn). (12) 
or 
r,=(exp(—Bn)), (13) 


where we used the expansion 


2 (1) k 
exp( pn)= ( = ; (14) 


Thus, the value of 7, is given by 


r= exp(—Bn)a,. (15) 
n=1 


Applying the approximate formula, 


> fn)= i "f(x)dx, (16) 

n=1 0 

we have: 

r= | exp(—Bxda(xax, 7) 
0 


or using the mean value theorem: 


rs=exp[— Bx()]. (18) 


We also used the integral form of the normalization condi- 
tion in Eq. (5): 
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[coax (19) 


0 


Note that we did not specify any specific laws of scattering 
probability in the derivation of Eq. (10). Thus, Eq. (18) can 
be applied in a much broader context than just scattering of 
photons by particles in a turbid medium. Comparison of 
Eqs. (12) and (18) shows us that ¥-(n) as 6—0. How- 
ever, generally speaking, x#(n). This is due to the fact that 
differences oj=(n!)—(n)! (j=2,3...°) are nonzero in 
general. 

The problem we face now is the determination of the 
parameter x. For this we use the well-known asymptotic 
result of the radiative transfer theory!”: 


B 1/2 
r,=1-4|.——_| , 20 
which is valid as 6—0. Here 
1 (r 
=; | cos Op(@)sin 6d 6, (21) 
0 


and p(@) is the phase function, which describes the prob- 
ability of photon scattering in a particular direction, speci- 
fied by the scattering angle 0. It is normalized as follows: 


1 T 
zl p(6)sin 0d0=1. (22) 


Clearly, for equal scattering probabilities: p(@)=1 and g 
=0. For phase functions highly elongated in the forward 
direction, which is the case for particles in atmosphere and 
ocean, g is close to 12 

Equation (18) takes the following form as B—0: 


ra l= £X. (23) 


So, comparing Eqs. (20) and (23) we have 


4 
x= ke? (24) 
where 

k=[3(1-g)p]'? (25) 


is the diffusion factor of the radiative transfer theory. Note 
that ii the diffusion approximation we have, instead of Eq. 
(25)": 


ka=[3(1 — gw) B]'”. 


The difference between k and k, is small in this case, 
which we consider here (wg— 1). 

Equation (24) shows that the average number of scatter- 
ing events in a semiinfinite layer is proportional to k~! as 
8—0. In particular, we have: x—-% as B—0 or g— 1. 

Finally, we have from Eqs. (18) and (24): 


s) 06 


r= ex = 


Equation (26) presents the total reflectance r, of a semiin- 
finite weakly absorbing medium as a function of two pa- 
rameters. They are the probability of photon absorption £ 
and the diffusion factor k. 

The value of k determines the decrease of photon flux in 
deep layers of absorbing semiinfinite layers. Namely, it fol- 
lows for the flux F at the depth z”: 


F=F o exe = Z| (27) 
La 

where Fy is the incident light flux and 

La=(kOex) | (28) 


is the so-called diffusion length. Note that Eq. (27) cannot 
be applied for nonabsorbing media.” The value of ox, gives 
us the extinction coefficient. Introducing the absorption 
length Labs, which is inversely proportional to the absorp- 


tion coefficient O41, 


Lass = ( Tabs) = > (29) 


we have from Eq. (26): 


< 


r.=exn| =7 (30) 


abs 


Thus, the total reflectance is determined by the ratio of the 
diffusion length L, to the absorption length Lavs: 


g= (31) 


abs 


For small values of & the value of r, is close to 1 and we 
have approximately: 


l-r, 


a oa (32) 


This value could be easily found experimentally for a given 
turbid medium. 

Here we come to the important conclusion. Namely, Eq. 
(26) can be applied to a much broader class of scattering 
processes than those described by the radiative transfer 
equation. To make things more clear, let us consider 
foamed media (e.g., whitecaps over ocean). The notion of 
the phase function can be hardly introduced for such a 
densely packed system of air bubbles. However, both the 
value Lps and Ly can be defined and in principle measured 
[see Eqs. (27) and (30)]. The application of Eq. (26) to the 
interpretation of light reflectance measurements for snow 
and foam samples has shown a high accuracy of the theo- 
retical model proposed.® This is mostly due to the fact that 
the photon free path length both in snow and foam is much 
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linear approximation 
exponential approximation 
=== exact computations 


Fig. 1 The dependence r,(y), calculated with Eqs. (20) (linear ap- 
proximation) and (33) (exponential approximation) and the exact ra- 
diative transfer code.’ 


larger than the wavelength. Then the conventional radiative 
transfer equation applies (with modified parameters of an 
elementary volume). 

Let us rewrite Eq. (26) in the more explicit form: 


1- wo 
3(1=8) 


This equation is often used in the applied radiative transfer 
theory.’ It can be applied to a much wider class of absorb- 
ing media than the simple linear approximation of Eq. (20). 
We notice that the combination of local optical characteris- 
tics, given by 


-4 1- wo 
Y3 (=g) 


is of paramount importance in the radiative transfer. It de- 
termines the total reflectance from a weakly absorbing 
semiinfinite layer. 

The accuracy of Eqs. (20) and (33) as compared to cal- 
culations with the numerical solution of the radiative trans- 
fer equation for water clouds with different radii of par- 
ticles is presented in Fig. 1. Specifically, we used the case 
of water droplets size distribution, given by f(a) 
= Na exp(—9a/a,;), where the effective radius ap was as- 
sumed to be equal to 4, 6, and 16 um. Here N 
=9/[r(a}] and T (7) is the gamma function. Calcula- 
tions were performed at a wavelength of 1.55 jm. The 
refractive index of droplets was assumed to be equal to 
1.3109— yi, where y was varied in the range 0.00001 to 
0.001. The values of g and wọ were found using the Mie 
theory. Exact results for all y and a. used in calculations 
lay on the same curve plotted against y. This is a well- 
known feature of radiative transfer in semiinfinite media, 
having the same value of y. 

We see that the accuracy of Eq. (33) is much better than 
that of an exact asymptotic result as in Eq. (20). Also it 


1/2 
r =ex| 4 | ; (33) 


1/2 
(34) 
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follows that Eq. (33) can be used at least up to y= 1/2 
regardless actual values of g and wọ. It means that the 
smaller the value of g, the larger the applicability of Eq. 
(33) in terms of the probability of photon absorption £. In 
particular, we obtain from Eq. (34) at y = 1/2: 


3 
Pin=zg (78). (35) 


Here ın is the value of 6, correspondent to y= 1/2 at a 
fixed g. We see that the value of B1n=3/64 (or approxi- 
mately 1.6%) at g=0. It is much smaller for values of g, 
which are close to 1. For instance, 64» is only 0.07% at 
g=0.85, which is typical for water clouds. Note that Eq. 
(33) can be applied for larger y (e.g., y e [1/2;1], see Fig. 
1), if a high accuracy is not a primary requirement (e.g., for 
rapid estimations). 

Our derivation underlines assumptions, which lead to 
the approximation of Eq. (33). Note that Eq. (18) follows 
directly from Eq. (13), assuming that o;=(n') —(n 0 Hor 
all j and x=(n). Then, 


(exp(— Bn))=exp(— B(n)). (36) 


This is in fact a major assumption used in the derivation of 
Eq. (33). 


3 Plane Albedo and the Reflection Function 


The same line of reasoning can be applied to the derivation 
of the solution for the plane albedo r,( 00), defined as 


1 27 a/2 
r(0)= 52 | av| sin 2 8d BR( Ty, O, Y). 
27 Jo 0 


Again, we have for an absolutely white Lambertian surface: 
R=1 and, therefore, =l for all incidence angles Jo. 

The MacLauren expansion of r,(0o) with respect to wo 
has the following form: 


ioe) 


ro) = 2 b,( 9) or, (37) 


where unknown coefficients b,, depend now on the inci- 
dence angle 3). Equation (20) takes the following form for 
the plane albedo”: 


r,(99)=1—uo(Bo)y, (38) 


where y is given by Eq. (34) and uo(%p) is a so-called 
escape function’ of the radiative transfer theory. It can be 
approximated by the following formula outside the region 
of the grazing incidence angles 99!: 


3 
uo( Bo) = 7 (1 +2 cos Do). (39) 


Equation (39) describes also the angular variation of light 
transmitted by overcast clouds. Then ®g in Eq. (39) should 
be substituted by 3.7 
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Following the same steps as before, it follows that, in- 
stead of Eq. (18): 


r)( 9) =exp[ — Bx(9o)], (40) 
where 

4u,(d 
reao) (41) 


Also we can write: 


rp( 8) =expl — yuo( Fo) ]. (42) 


It follows from Eqs. (39) and (41) that the average number 
of scattering events x(0 9) as B—O is larger for smaller 
angles o. This effect follows for physical grounds as well. 
Indeed, it is harder for photons to escape the medium if 
they are injected along the normal to the scattering me- 
dium. We see, therefore, that the contribution of multiple 
scattering decreases with Jo. 

The final task is to obtain the similar approximate equa- 
tion for the reflection function of a semiinfinite medium 
R..( 9,9, ).” We use the following MacLauren expansion 
in this case: 


oo 


R99, 9, W= Dy c(d, 0, Po. (43) 


n= 


R (Do, Ò, p) is not necessarily equal to one at wọ=1. So, 
we have: 


(os) 


REC 90 D= caldo F 1, (44) 


n= 


The subscript 0 means that the quantity RÈ? corresponds to 
the case of zero absorption [R2=R.,()=1)]. Let us con- 
sider the ratio 


R..( o, O, y) 
ï= a (45) 
R( Do ? Ò, 0) 
It can be represented by the following expansion: 
Y= Di dul 80,9. Ho, (46) 
where 
Cal Vo > 0, p) 
d,(®o, È, Y) = se l (47) 
° v =e nos. Y) 
It follows that 
d (ðo, ®, Y)=1, (48) 
l 


n= 


and numbers d,, can be interpreted in terms of probabilities. 


Repeating the previous steps, we arrive at the following 
formula: 


Y=exp[— BX(, 3, )]. (49) 


Again we use the exact asymptotic radiative transfer theory 
to get the average scatterings number x( 0, 0, Y). Namely, 
it follows that”: 


Y=1-h(9), Wy, (50) 
where 

0 ð 
i. 8. a, 61) 


On the other hand, we have from Eq. (49): 
Y=1-Bx(9, È, y). (52) 


Therefore, comparing Eqs. (50) and (52), we arrive to the 
following important equation: 


4h( do, O, 
naapa A (53) 


Finally, it follows [see Eqs. (45) and (49)] that: 
RC 8,9, P) =RA( 80, 9, Wexpl—yh(G, 9H], (54) 


where we used the equality: y=46/k. 

Formulas (33), (40), and (54) have been known for a 
long time.” However, the new way of deriving them can 
facilitate applications. Also we note that the function h, 
given by Eq. (51), has a simple physical sense. It is equal to 
the ratio X( 99, 0, w)/x [see Eqs. (53) and (24)]. Therefore, 
the function h gives the ratio of the average number of 
scatterings for a given incidence, observation, and azimuth 
angles to that for the case of total reflectance. 

Results presented previously can be generalized for the 
case of polarized light. Then the Stokes vector parameter 
with components J, Q, U, and V> should be used instead 
of diffused light intensity. In particular, let us consider the 
case of the normal illumination of a semiinfinite plane- 
parallel turbid layer. Then the degree of polarization of re- 
flected light is given by: 


P.=- — (55) 


and components U and V vanish. The positive polarization 
means that oscillations of the electric vector are preferen- 
tially in the plane perpendicular to the plane containing the 
observing direction and the normal to the layer. 8 

It is known that the value of J is much more sensitive to 
the change of the absorption in the medium than the param- 
eter Q. To simplify, we assume that Q in Eq. (55) coincides 
with its value for a nonabsorbing medium. 

To study the validity range of this assumption, we have 
calculated values of Q and J using the doubling method of 
the vector radiative transfer equation solution, assuming 


Optical Engineering, Vol. 42 No. 7, July 2003 5 


Kokhanovsky: Optical properties of semiinfinite . . . 


that the optical thickness 7 is equal to 500, which is a close 
representation of a semiinfinite medium. The wavelength 
was equal to 0.65 um. The phase function was kept con- 
stant and equal to that of water droplets with the same 
particle size distribution as in Fig. 1 (at a,s=6 wm). The 
value of the asymmetry parameter of g was, therefore, con- 
stant and equal to 0.85. The single scattering albedo was 
changed in the range 0.95 to 1.0. This corresponds to the 
change of the probability of photon absorption 6 from 0 
until 5%. Then the parameter y is changed in the range 
[0,4/3]. 

It should be stressed that such small variations of the 
single scattering albedo produce large variations of the 
value of Z and the reflection function”: R = al/ HoE. Here E 
is the incident light flux density. We found, in particular, 
that 7 changes more than four times at the nadir observa- 
tion, and a 60-deg observation angle for changes of albedo 
in the range 0.95 to 1.0 (a semiinfinite medium case). The 
value of Q changes only by 13% for the same conditions. 
The variation of Q is below 5% for photon absorption prob- 
abilities below 2%. 

Therefore, assuming that values of Q for nonabsorbing 
and weakly absorbing media coincide, we obtain, using 
Eqs. (54) and (55): 


P..( 9) = PL 9)explyh(0,8,0)], (56) 
where P? =-— R21 /R}, and we used equalities 
I=Rl0, Q=Rozl0, (57) 


which are valid in the case of incident unpolarized light 
having the intensity J). Here R2; is the element 21 of the 
reflection matrix of a turbid medium. We neglect the de- 
pendence of R2; on y and assume that Razi RY aj 

Equation (56) shows that the increase in absorption 
leads to an increase in light polarization, which is a well- 
established result.!° Moreover, the ratio P?IP, coincides 
with Y [see Eqs. (49) and (56)]. 

Errors in Eqs. (54) and (56) are given in Fig. 2, assum- 
ing the same particle size distribution of water droplets as 
in Fig. 1 at ay=6 wm, \=0.65 um, By=0 deg, and 3 
=60 deg. The single scattering albedo was changed in the 
range 0.95 to 1.0. We found that the error is smaller than 
10% at wg>0.98 (R..>0.32) for Pa [see Eq. (56)]. It is 
smaller than 8% at wy>0.95 (R..>0.18) for the value of 
R [see Eq. (54)]. 

Therefore, approximations presented here can be used to 
estimate the degree of polarization and intensity of light 
reflected from weakly absorbing (@ g— 1) semiinfinite (7 
—o) media, avoiding the numerical solution of four 
coupled integrodifferential radiative transfer equations,® 
which could be of importance for optical engineering ap- 
plications. 


4 Conclusions 


We derive a number of approximate equations, valid for the 
case of semiinfinite (r—œ) weakly absorbing (wọ—> 1) 
turbid media. They can be used in a wide range of applica- 
tions, especially for the development of the spectroscopy 
and polarimetry of optically thick, weakly absorbing turbid 
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reflection function 


error, % 


0.9 0.96 0.97 0.98 0.99 1.00 
single scattering albedo 


Fig. 2 Errors of Eqs. (54) (solid line) and (56) (dashed line). Exact 
calculations were performed using the vector radiative transfer code 
for the same particle size distribution of water droplets as in Fig. 1 at 
āe 6 wm, \=0.65 um, 3p>=0 deg, and d=60 deg. 


media. One important limitation is the infinite thickness of 
the layer assumed, which is not always the case (e.g., for 
clouds). However, it should be pointed out that there are 
simple equations, which connect reflection and transmis- 
sion characteristics of turbid layers with optical thicknesses 
larger than 5 with functions derived here [e.g., see Eq. 
(54)].17®1!-!3 This extends the results presented here to a 
broader range of turbid media, including clouds, leaves, 
paper, and photolayers. 

Equations (33), (42), (54), and (56) can be used for the 
determination of the spectral dependence of a single scat- 
tering albedo of the scattering medium in question, which 
allows us in principle to find the spectral behavior of the 
imaginary part of the refractive index of scatterers. This can 
be used for spectroscopic purposes and the chemical com- 
position of particles identification.” 

It follows from Eqs. (56) that darker surfaces yield 
larger values of the degree of polarization. This is con- 
firmed by experiments in Refs. 14 and 15. 
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